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SUMMARY 

Jet-A spray combustion has been evaluated in gas turbine combustion with 
the use of propane chemical kinetics as the first approximation for the chemi- 
cal reactions. In this work, the numerical solutions are obtained by using 
the KIVA-II computer code. The KIVA-II code is the most developed of the 
available multidimensional combustion computer programs for application of the 
in-cylinder combustion dynamics of internal combustion engines. The released 
version of KIVA-II assumes that 12 chemical species are present; the code uses 
an Arrhenius kinet ic-control led combustion model governed by a four-step global 
chemical reaction and six equilibrium reactions. Our efforts involve the addi- 
tion of Jet-A thermophysical properties and the implementation of detailed 
reaction mechanisms for propane oxidation. Three different detailed reaction 
mechanism models are considered. The first model consists of 131 reactions 
and 45. species. This is considered as the full mechanism which is developed 
through the study of chemical kinetics of propane combustion in an enclosed 
chamber. The full mechanism is evaluated by comparing calculated ignition 
delay times with available shock tube data. However, these detailed reactions 
occupy too much computer memory and CPU time for the computation. Therefore, 
it only serves as a benchmark case by which to evaluate other simplified 
mode 1 s . 

Two possible simplified models have been tested in the existing computer 
code KIVA-II for the same conditions as used with the full mechanism. One 
model is obtained through a sensitivity analysis using LSENS , the general 
kinetics and sensitivity analysis program code of D.A. Bittker and 
K. Radhakrishnan. This model consists of 45 chemical reactions and 27 species. 
The other model is based on the work published by C.K. Westbrook and 
F.L. Dryer. This simplified model consists of 5 chemical reactions and 12 
species. The numerical results indicate that the variation of the maximum 
flame temperature is within 20 percent as compared with that of the full mech- 
anism of 131 reactions. The chemical compositions of major components such as 
CqHq, H-?0, Oo, CO 2 , and N 2 are of the same order of magnitude. However, the 
concentrations of pollutants are quite different. Details are presented later 
in this paper. 


INTRODUCTION 


The goal of gas turbine combustion system design and development is to 
meet many mutually conflicting trade-offs between high combustion efficiency 
over a wide operating envelope and low N0 X emissions: low smoke, low lean 

flame stability limits, and good starting characteristics; and low combustion 
system pressure loss, low pattern factor, and high structural durability. 

Recent advances in high-speed and large-memory computers, numerical algorithms, 
grid generation, and combustion modeling have a powerful positive influence on 
future design capability. To accurately model combustion system subcomponents, 
significant improvements in physical modeling and numerical techniques are 
required. 

In this work, a detailed chemical kinetics mechanism of a simple hydrocar- 
bon, such as propane, is used to predict combustion flows and emissions of 
typical combustion systems representative of gas turbine combustors using 
multidimensional codes. Previous multidimensional computations of combustion 
and flow characteristics of gas turbine combustor sector or subcomponents used 
global or quasi-global chemical kinetics mechanisms. These combustion models 
may be able to predict qualitatively the average value of heat release; how- 
ever, they are unable to predict local chemical equilibrium and detailed com- 
bustion species profiles, including N0 X emission. In the current work, the 
fuel used is Jet-A, whose physical properties are known and are given later in 
this paper. However, the chemical kinetics for Jet-A combustion with air are 
not presently available. It is necessary to use the propane chemical kinetics 
for the simulation of Jet-A combustion as a first step approximation. However, 
even the chemical kinetics for propane are not well understood. In this work, 
the propane oxidation mechanisms consisting of 131 reactions (ref. 1) are used 
to perform benchmark computations and serve to evaluate the other two reduced 
mechanisms . 

The numerical simulation is done by using the code KIVA-II, which was 
originally designed for the simulation of turbulent reactive flows in piston 
engines. The code is modified for testing propane chemical kinetics in two- 
dimensional gas turbine combustion. The purpose of using the KIVA-II code is 
that it is considered to be very well developed and represents present state- 
of-the-art. A description of KIVA-II and its modifications will follow. 


DESCRIPTION OF COMPUTER CODE KIVA-II 

The computer code KIVA-II was originally designed to simulate numerically 
transient two- and three-dimensional reactive fluid flows with sprays in piston 
engines. The dynamics and evaporation of spray droplets, spark ignition, the 
chemical kinetics of combustion, turbulence, and the boundary effects of the 
flow near the walls are considered in detail. Turbulence is modeled by using 
two-equation k-e model equations and an implicit-continuous Eulerian (ICE) 
technique is used for the flow solver. A stochastic particle method is used 
to calculate the evaporation rate of liquid sprays. The effects of droplet 
collisions and aerodynamic breakups are included in the computations. Octane 
was the fuel originally used in the combustion. Properties of several other 
hydrocarbon fuels are provided in the subroutine fuel, so that the combustion 
of other fuels can be simulated. However. Jet-A properties are not included. 
Detailed descriptions of the equations and the numerical method of solution 
are discussed in reference 4. 
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MODIFICATION OF KIVA-II 
Jet-A Fuel Properties 

In order to use Jet-A as the fuel for combustion, the properties of Jet-A 
were implemented in the subroutine fuel. The information includes the 
fol lowing. 


Chemical symbol, Cj 2 H 23 

Molecular weight, 176.315 

Critical temperature, 668 K 

Heat of formation (0 K), 72.53 Kcal/raole 

Surface tension (350 K), 18.85 dynes/cm 

The enthalpies, latent heat, vapor pressure, and liquid viscosities of 
Jet-A fuel are given in appendix A from a computer printout. 


Modification of Subroutine "RINPUT" 


Because the number of species appearing in the combustion of propane is 
different from that of the combustion of octane, the enthalpy data of species 
must be modified in the subroutine RINPUT. The numbers of species for 131, 
45, and 5 chemical reactions are 40, 27, and 12, respectively. The enthalpy 
data in the subroutine RINPUT are arranged according to the order of species 
appearing in the reactions. All the data are listed in the temperature range 
from 0 to 5000 K in increments of 100 K. A list of the 40 species is given 
below. 


C 3 H 8 , 0 2 , n 2 , C0 2 , h 2 o, H, H 2 . 0 , N 
OH, CO, NO, CH, CH 2 , CH 3 , CH 2 0, CH 3 0 

c 2 h 3 , c 2 h 5 , c 2 h 6 , c 3 h 7 hco, ho 2 , no 2 , ch 4 , c 2 h 2 
C 9 H 4 , CN, c 2 h, c 2 ho, hcn, hnco, hno, hno 2 
HN0 3 , H 2 0 2 , NC0 , NH, NH 2 , N 2 0 

A list of enthalpies for the 40 species appears in appendix B. 


Modification of Input Data File 

The major change is in the data provided for the chemical kinetics. As 
an example, the input file for the simplified chemical mechanism is shown in 
appendix C. 


JUSTIFICATION OF THREE MECHANISMS FOR THE COMBUSTION OF PROPANE WITH AIR 

The full mechanism, consisting of 131 chemical reactions, is considered 
to be the standard for this evaluation. This standard is established through 
two extensive studies. In the first study (ref. 1 ) the numerical predictions 
of pollutants NO and CO were compared with experimental measurements. They 
were found to be in good agreement. In the second study (ref. 2 ) the number 
of reaction equations was already 131 . The full mechanism was further sup- 
ported by the ignition delay time checked with the shock tube data. 


3 



The mechanism with 45 reactions has been considered to be the first possi- 
ble reduced mechanism for saving computer memory capacity and CPU time. The 
45 reactions were determined through sensitivity analysis by mixing propane 
and air in a container with a definite volume. The results of temperature, 
pressure, and species concentrations versus time from 45 reactions were com- 
pared with those from the full mechanism. They were in fairly good agreement 
(ref. 2). 

The combustion of propane in the simplified mechanism consists of five 
reactions. Two of the equations are obtained from Westbrook and Dryer 
(ref. 3), and the other three are from the original input file of KIVA-II 
describing the extended N0 X Zeldovich mechanism. In addition, the six equilib- 
rium reactions described in reference 5 are also used with this five-kinetics- 
reaction mechanism. The five reactions, together with the reaction constants, 
are given below. 

1. C 3 H q + (3/2 + 8/4) 0 2 - 3C0 + 4 H 2 0 

k f - 10 12 Exp(-30/RT) [C 3 H 8 J° 1 [t^] 1,65 
K b - 0 

2. CO + 1/2 0 2 = C0 2 

K f = 10 14 - 6 Exp(-40/RT) [CO] 1 [H 2 0]° * 5 [0 2 ]° 25 
K b = 5xl0 8 Exp(-40/RT) [CO^ 1 

3. 0 2 + 2N 2 =. 2N + 2N0 

K f = 1.5587xl0 14 Exp(-6.7627xl0 4 /T) [0 2 ] [N 2 ] 2 
K b = 7.500xl0 12 Exp ( 0 ) [N] 2 [NO] 2 

4. 20 2 + N 2 = 20 + 2 NO 

K f = 2.6484xl0 10 Exp(-5 . 9418xl0 4 /T [0 2 ] 2 [N £ ] 

K b = 1.6xl0 9 Exp(-1.9678xl0 4 /T) [0] 2 [NO] 2 

5. N 2 + 20H = 2H + 2 NO 

K f = 2.123xl0 14 Exp(-5.7020xl0 4 /T) [N] 2 [OH] 2 
K b - 0 
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RESULTS AND DISCUSSION 


Comparisons With Data of Direct-Injection Strat i f ied-Charge 


Piston Engine Combustion 

The three mechanisms of propane chemical kinetics were tested first in 
KIVA-II with computations of direct-injection strat if ied-charge (DISC) combus- 
tion in gasoline- fueled piston engines. The specifications and test conditions 
for the piston-cylinder configuration are listed in table I. For these calcu 
lations , the numerical grid consists of 20 grid points in the axial direction 
and 22 grid points in the radial direction. A Bessel function velocity profile 
(ref. 4) was used to simulate the swirl velocity profile of the incoming air. 
Standard boundary conditions were used in this two-dimensional axisymmetric 
engine model and the details are described in reference 4. 

The three mechanisms of propane chemical kinetics were evaluated by com- 
parisons of data obtained from the numerical simulations of the DISC piston 
engine. In all test cases, the simulations were started at the beginning ot 
the compression stroke and ended at the crankangle of 65° ATDC. For the base- 
line case, case 1, a quasi-global mechanism of isoctane oxidation was used 
together with the extended N0 X Zeldovich mechanism and six equilibrium reac- 
tions (ref. 5). The propane chemical kinetics mechanism for case 2 consists 
of five kinetic reactions and six equilibrium reactions (these six equilibrium 
reactions are the same reactions considered in case 1). In addition, the 
45-reaction mechanism and the 131-reaction mechanism of propane oxidation were 
also examined (cases 3 and 4). Table I summarizes the test conditions tor 
cases 1 and 2. For cases 3 and 4, the starting crankangle of ignition was set 
at 27° BTDC for an ignition duration of 9.6 degrees. In all cases studied, 
liquid isoctane was injected starting from 52° BTDC for an injection duration 
of 12.672 degrees in an amount that corresponded to an overall equivalence 
ratio of about 1.0. 


Figures 1 and 2 show the computed mean gas velocity, fuel spray, tempera 
ture profile, fuel vapor, and CO profiles at a crankangle of about 29° and 13 
BTDC for cases 1 and 2. Comparisons of these results indicate that the tive 
reaction mechanism of propane oxidation (case 2) predicts a similar combustion 
rate and slightly lower maximum flame temperature compared with the results 
predicted by the four-reaction mechanism of isoctane oxidation (case 2). These 
results also indicate that the flame temperatures are lower than the corre- 
sponding local equilibrium temperature. 


Figures 3 and 4 show the computed mean gas velocity, fuel spray, tempera 
ture, and fuel vapor profile near TDC for cases 3 and 4. It may be seen that 
the flame is not well propagated and diffusive throughout the engine chamber. 
The temperature plots of figures 3 and 4 reveal that the predicted maximum 
flame temperature concentrates only within a small region downstream ot the 
injector tip. This observed low combustion efficiency is probably due to poor 
fuel-air mixing, ignition, and combustion rate as a result of non-optimal 
engine design and test conditions. The maximum flame temperatures observed 
for cases 3 and 4 are lower than the corresponding equi librium. values. The 
observed flame temperature from the mechanism of 45 reactions is lower than 
that of the full mechanism. This indicates that in the 45-reaction mechanism 
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more endothermic reactions and fewer exothermic reactions are considered on a 
relative basis. 


Comparisons With Data of Model Airblast Fuel Nozzle/Reactive Flow 

The model airblast-type fuel nozzle used in this study is shown in fig- 
ure 4. It is a typical fuel nozzle used in gas turbine combustors. It embod- 
ies inner and outer airstreams which help to atomize and distribute fuel 
sheared from a thin annular sheet. The model fuel nozzle injected fuel bv 
means of an annular fuel passage with the nozzle tip exit area of 0.03 cnr . 

The fuel exit area is defined by the tip of the core inner air pipe and the 
position of the wedge-shaped fuel-filming lip of the fuel flow passage. No 
swirl is imparted to the fuel flow. This fuel nozzle produces a hollow-cone 
fuel spray with a nominal included angle of 75 degrees. The nonswirl inner air 
is supplied to the core air pipe with exit an diameter of 1.0 cm. The outer 
airflow is swirled by a vane assembly. Different vane assembly designs were 
evaluated. The end cap of the outer air shroud imparts a radial-inflow compo- 
nent to this airflow; the end cap exit area is 3.24 cm^ . 

The flowfield downstream of the fuel nozzle is axisymmetric without con- 
finement. For these calculations, the numerical grid consists of 60 grid 
points in the axial direction and 52 grid points in the radial direction. 
Boundary conditions are applied to the inlet, outlet, solid walls of the fuel 
nozzle, and centerline of the unconfined flow downstream of the fuel nozzle. 

A brief description of these boundary conditions are as follows. In this work, 
the effect of outer air swirler designs on fuel atomization and spray patterns 
is examined. For the baseline case, the outer air inlet flow is assumed to 
have a zero mean radial velocity and is characterized by a uniform axial veloc- 
ity and pressure and an azimuthal velocity profile of a solid body rotation. 

The inner air is assumed to have zero mean radial and azimuthal velocity, and 
enters the core air pipe with uniform axial velocity and pressure. The temper- 
ature of the inner and outer air is assumed uniform and equal to 1367 K. The 
profiles of turbulent kinetic energy and turbulent length scale are taken from 
reference 6. At the solid walls the velocity component normal to the wall is 
set to zero while the law of the wall is applied to the velocity component 
parallel to the wall. The normal gradient of turbulent kinetic energy near the 
wall is set to zero and the dissipation rate of turbulent energy is calculated 
by using an equilibrium boundary layer. The turbulent length scale is assumed 
to vary linearly with distance from the solid walls. At the unconfined flow 
centerline both the radial and azimuthal velocity components are zero, while 
the radial derivatives of the other flow quantities are zero for reasons of 
symmetry. The model airblast fuel nozzle is sized for nominal inner air, outer 
air, and fuel flow rates of 1.3 gm/sec, respectively. A Rosin-Rammler droplet 
distribution function is used. The calculations are performed with a time step 
corresponding to 5xl0 _ ® sec, and 3xl(P droplet parcels are injected at each 
time step. 


CONCLUSIONS 

Throughout this work it can be seen that the numerical simulation of 
octane combustion in KIVA-II is used as a guide for propane combustion because 
accurate experimental data for octane combustion are well known and available 
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for direct comparison for cascade validation. On the other hand, the adiabatic 
flame temperature of propane combustion is another checkpoint to be considered. 
When the air temperature is 1033 K, the propane temperature 311 K and the pres- 
sure 1 atm, the adiabatic flame temperature is 2560 K. Because of the vapori- 
zation of liquid fuel, either Jet-A or octane, and the heat transfer to the 
surroundings, the flame temperature is lower than the adiabatic flame 
temperature . 

All simulations are carried out on the CRAY X-MP computer at the NASA 
Lewis Research Center. Because the chemical kinetics are only a part of the 
whole simulation, the changes in the CPU time and memory space which can be 
directly attributed to the kinetic mechanisms for three different mechanisms 
cannot be clearly accounted for at the present time. Consequently, no 
definite numbers are reported. 

In the prediction of pollutants N0 X and CO, it is observed that the full 
mechanism of 131 reactions is more reliable than the other two mechanisms. 

Under the same conditions stated above, a comparison between the combus- 
tion of octane and propane is made. The profiles of velocity and temperature 
are presented in figures 5(a) and (b) . These simulations are performed with 
the understanding that, for the first 1000 cycles of calculations, the combus- 
tor is filled only with the airflow, the fuels are injected after 1000 cycles 
and then ignited. The ignitor is located on the fuel path. The amount of fuel 
is nearly stoichiometric. The chemical kinetics for the combustion of octane 
consist of 4 reactions, and for propane 5 reactions. The mechanisms of 131 and 
45 reactions are to be investigated in the near future. 
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APPENDIX A 

JET- A FUEL PROPERTIES 


C + + + 

C ++ + THE ENTHALPY DATA BELOW ARE COMBINED FROM JET-A(L) 

C *+« JET-A (CI2H23) IN THE TEMPERATURE OF C — 500K AND JET-A(G) IN THE 
C ++♦- TEMPERATURE OF 600 — 5000K 

data ( h k ja ( n ) , n« 1 * 51 > / 0.00, 7. 362, 1 5. 342, 22. 692, 32. 002, 43 272. 

1 34. 21. 66. 38, 79. 61, 93. 77, 108. 71, 124. 30, 140. 44, 157. 00. 174. 13, 

2 191. 61, 209. 42, 227. 53. 24 5. 90, 264. 51, 283. 31, 302. 29, 321. 42, 340. 68, 

3 360. 06, 379. 52, 399. 07, 418. 69, 438. 37, 450. 11, 477. 89, 497. 72, 517. 59, 

4 537. 30, 557. 44. 577. 42. 597. 44, 617. 49, 637. 50, 657. 70, 677. 85, 698. 05. 

5 710. 27, 738. 52, 750. BO. 779. 10, 797. 41, 019. 73, 040. 04, 860. 34, 880 62/ 

c ^ 

c ++*» the liquid latent heats, in ergs/gm, at intervals t=100(n-l): 
c + + + values between 300k and upper limit taken from n. b. vargaftik 
c + + > '’tables on the thermophysical properties of liquids and gases" 
c +++ values are linearly extrapolated back to t=0 k. 

c + + + 

c + + + 

c +++■ Jet-A latent heat values in range 0-500k : 
c +++ 

data (hlatja(n), n=*l, 51 ) /3. 792e9, 3. 598e9, 3. 405e9, 3. 208e9, 

1 2. 842e9, 2. 425e9, 45*0. 0/ 

c +++ 

c ++.»- 

c + + the liquid vapor pressures, in dynes/sq. cm. , at intervals 

c n + t-lO(n-i). also taken from vargaftik. note that each table 

c 4- + ■*■ must contain values beyond tcrit of the fuel to protect against 

c +++ over-reaching in subroutine evap. also note that pvap is 

c +++ dimensioned at lvap words in common /lek/ in comdeck comdi 
c •**■**• lvap must be adequate for all fuels in the library, 
c +++ for heptane, dodecane, tridecane and 

c +++ hexadecane. the vapor pressure is scaled with temperature 
c + + + using eqn (see reid, prop, of liq. and gas. eqn 6—2. 4,2. 5) 
c +++ 1 n ( pv ) -1 n ( pvO )»( I n < pvl ) -In CpvO) )(< t-tO )/( 1 1 -tO >) *t 1 /t 

c +++ w / tl=tcrit, t0=520k (dodecane ), t0=510k (tridecane) 

c f + + 

c + + + Jet-A vapor pr#»sur« avenj lOk in ranjje 10-690k ( ter i t=A68k ) : 

C + + -*■ 

data (pvapja(n), n=l, 69) /. OOE+OO, . OOE+OO, . OOE+OO, . OOE+OO, . OOE+OO, 

1 OOE **00, . OOE+OO, . OOE+OO, . OOE+OO, . 29E-09, . 18E-07, . 57E-06, . 1 IE-04, 

2 . 13C-03, . 1 1 E— 02, . 76E-02, . 4 1E-01, . 18E + 00, . 6QE +00, . 23E*-01, . 67E+01, 

3 . I8E+02, . 44E+02, . 1 0E+03. . 22E+03, . 44E+03, . 04E+O3, . 15E>04, . 27E+04, 

4 . 4 5E »-04 , . 7 4 E+04 , . 12E+05. . 10E+O5, . 27E+05, . 40E+05, . 57E+05, . 80E + 05, 

5 I 1 E + 06 » . 15E+06, . 20E+06. . 27E+06, . 35E+06, . 45E+06, 57E+06. . 72E+06, 

6 . 90E *06, . 1 1E + 07, . 14E + 07, 16E+07, . 20E + 07, . 24E+07, . 2BE+07, . 33G+07. 

7 . 39C+07, . 45E +07, . 53E + 07. . 61E+07, . 70E+07, . 00E+07, . 91E+07. . 10E+08. 

8 . 1 2E+08, . 13E+08, . 1 5E+0B, . 16E+0B, . 18E+08, . 2OE+O0, . 22E+08, . 24E + 0B/ 

c + + + 

c + + + the liquid viscosities, in gm/(cm sec). as in the vapor pressures 

c •+ + + above, the intervals are t = 10(n-l), and the tables contain values 

c +4-+ beyond tcrit of the fuel. note the visual e x trap o 1 a t i on required 
c + + + for dodecane, tridecane, tetradecane and hexadecane, for' lack of 
c + + + tabular data. this needs improvement, esp. around tcrit. 
c + + + 

c + * + Jet-A liquid viscosity every 10k from 0-4I0k (tcrit=668k) data 
c + ++- obtained from HD of Aviation Fuel Properties 1984. 

data (vislqja(n), n — 1 # 4 2) /23*0. 14 02, 0. 0081 1,0. 05607. O. 03045, 

1 O 02804, O. 02163, O. 0 1722, O 01386, 0. 01161, 0. 00985, O. 00857, 0. 0074 5, 

2 0. 00657, 0. 00593, 0. 00529, 0 00473, 0. 00433, 0 00392, 0. 00360. O. 00336/ 


r 

c COO 0000000000000-00000 o oo 
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APPENDIX B 

ENTHALPIES OF ADDITIONAL SPECIES 


data 13).nM.51 >'-2. 

1 2. 12. 2.1 i. 5.5*. 4. 

2 1.15. 10. Or. 10.90, 11. 

1 17.54. 18.32, IV?8. 20. 

4 24.20. 2/. 21. 28.22. 29. 

5 55.54. 54.57. 17.59, 54. 


data Chkln. 14). nM . 51 )/-2. 

1 2.54. 5.*7. 4.45. 5. 

2 12.04. 15.24. 14.44, 15. 

5 25.44. 24.79. 24.15. 27. 

4 55.45. 37.00. 38.34. 39 

5 48.07. 49.47. 50.44, 52. 


.-1.581.-0,785. 0.015, 0.833. 1.679, 
4.44. 7,50. 8.60. 9.73. 10.84. 

14.98. 14,25. 19.54. 20.43. 22.14, 

24.82- 50 17, 31.53. 32.89. 34.26, 

41.14. 42 52. 43.91. 45.29. 46.68, 

53.46. 55.06. 54.46. 57.87. 59.28/ 


CM3 

da t * C Hk ( a . 1 5 1 . 

1 3.12. 4.31. 

2 15.47. 17.51. 

3 31.51. 33.33. 

4 44.24. 50.15. 

5 45.42. 67.34. 


nM. 51 >/“2. 477.-1 .482.-0.845, 0.017 , 0.97 
V 54. 6.44 . 8.25. 9.69. 1 1.17 . 12.70 

19.18. 20.44, 22.60, 24.34, 26.11. 27.19 

35.17. 37.02. 34.44. 40.74, 42.61, 44.49 

52^05. 53.95. 55.85. 57-76. 59.67. 41.54 
69.27, 71.20, 73.13. 75.06. 77.00, 74.94 


0.975. 2. 
2.70. 14 . 


1 2 99. 4 20 . 5.49, 6.87 . 4.31 , 9.42. 1 1.38. I?-™* *2 ■ 

2 14 29. 18. 00. 19.74. 21.50. 23.29. 25.09. 24.91. 2* 75, 30. 
» 44 34 34. 36 22, 34.11, 40.00. 41.90, 4^.30. 45.71, 47. 
I H'W. 53:37! 5 5.29. 57.22. 59.14. 6 1.07 , 6 3.00, 44 . 
5 64.47. 68.81, 70.75. 72-69. 74.64. 76.54, 74.54. 40.49. 82. 


lata (Kk(n.l7). 
3.48. 4.96. 

20.20, 22.38, 

40.94, 43.35, 

42.94, 65.46. 
85.42. 87.93. 


n*l.51 1 /- 2 . 431, -1.434, -0.834, 0 ■ 0 1 7 * 1 * ?? 5, , \ ' 

6.57, 8.24, 10.09, 1 1.99, 13.94. 15.96, JJ* 

24 40. 26.35. 29.14. 31.46. 33.80. 36- 16. -.8. 

45.77. 48.21, 50.64, 33.11. 55-52. 58.03. JO- 

67 94. 70.43, 72.92, 75-41. 77.91, 80.41, 82. 

90*44, 92.96. 95.48. 98.01,100.54,103.08.105. 


data (Hktn.U) 

1 3.61. 5.12 

2 20.39. 22.56 

3 41.03. 43.44 

4 65.01. 65.48 

5 85.41. 87.92 


. n*l. 51 >/-2. 458.-1 .66 5.-0. 460. 0.0 14, 1.046, 2. 
, 6.75. 8.48 . 10.29. 1 2.19. 14. U, U - 1« ■ J* ■ 
, 24.77. 27.01. 29.29. 31.59. 33. 9 t . 3 6 . 28 , 38. 
. 45.85. 48.28. 50.71. 53.16. 55.61. 58.07, 

47 96. 70 44. 72.93. 75.42, 77.91. 80.41. 82. 
90 43, 92.95. 95.47, 97.99,100.52,103.05.105. 


(hkfft,Zn.n*1.5l)'-?. 514.-1 719.-0.909. 0,019, 1.161. 2.515. 
13. 5.91. 7.85. 9.92. 12. in, 144ft. 14.77, 19.23, 21.76, 

.35, 26.99, 29.69. 32.43. 3S.22. 38 04. 40 .89, 43.77. 44.68. 

.61. 52.56. 55.52. 58.50. 61.49, 64.49, 47.50, 70.52, 73 55. 


DATA (HUM. 2* ) .11*1 . 51 >/-2 .073 

1 2.15. 2.89, 3.45. 4.42. 

2 9.35. 10.19. 11 .05. 11 .93, 

3 17.41. 18.36. 19.32, 20.29. 


* 1 . 51 >/-2 . 07 3, - 1 . 390, -0 .685. 0.013. 0,71 3. 1.422. 

3.6 5. 4.42. 5.22. 6.02, 4.83. 7 .46. 4 .49. 

11-05. 11.93, 12.81, 13.71, 14.42. 15-54, 16.47, 
19.32, 20.29. 21.27. 22.26, 23.26. 24.2 7. 25.29. 


4 26 31. 27.35. 28.39, 29.43. 30.49. 31.55. 32-42. 33.69. 34.77, 

5 35.85. 36.93. 38.02, 39.12, 40.21, 41.31. 42.41. 43.52. 44.62/ 


DATA ( HA( II, 29) ,M*l . 51 )/-2 . 356 

1 3.12. 4.24. 5.40. 6.40. 

2 14.62. 16.04. 17,53. 19.01. 

3 23 . 20. 29.76, 31 .32, 32.44 . 


MM . 51 )/-2 . 356 . -1 .600,-0 . 584. D.tll, 1.002. 2.04 1 , 
5.40, 6.40. 7.84. 9.12. 10.45. 11.81. 13.20. 

17.53. 19.01. 20.52. 22.03, 23.56. 25.10. 26 65. 

31.32, 32.54. 34.45, 36.01. 37.58, 39.14. 40.70, 

4 5.38, 46.93. 48.48. 50.03. 51-54. 53.12. 54.67. 

59.24, 60.42. 62. 35. 63.59, 65.42, 66.95. 68.48/ 


DATA ( HK ( 11 , 

30). MM ,51 )/ 

-2.636 

-1 . 888 

-1.004 

0.020 

1.149 

2 .359 

1 3.43. 

4 

. 97 . 6.36, 

7.81, 

9.31. 

10.84, 

12.44, 

14.07, 

15.72. 

2 17.41, 

l 9 

.12, 20.86, 

22.62, 

24.40, 

26.20, 

28.02, 

29.85, 

31 -69. 

3 33.54. 

35 

.41, 37.2S. 

39.lt. 

41 .04. 

42.94, 

44.85. 

46-74, 

48 . 64. 

4 50.55. 

52 

.47, 54.38, 

54 . 30 . 

58. 23. 

60.15. 

62.08. 

64.01, 

6 5.94, 

5 67.87, 

69 

.81, 71.75, 

73.69, 

75.64, 

77.58. 

79.53, 

8 L . 48. 

8 3 . 44/ 


DATA (HKII.31 >. MM. 51 >/-2. 207.-1 .513,-0.794. 0.016, 


13.74. 15.07. 


39.30. 40,77 
52.67, 54.16 


29 . 07 . 30 , 52 . 
42 , 25 . *3 73 . 


,-1.513,-0.794, 0.016, 0.915, 1.883, 

7.42. 8.44, 9.88. 11.14. 12,43. 

19.15. 20.54. 21.94. 23.35. 24.77. 

31.97, 33.43, 34.89, 36.36, 37.83. 

45.22, 44.70. 48,19. 49.68. 51-17. 

58.66. 60.17, 61.67. 63.18, 64.69/ 


DATA < HJUII. 32) . NM . 51 >/- 2 -620 , 


69.43. 71.38, 


-1 424,-0.980, 0.020. 1.167. 2.432. 

985. 11.48. 13.15. 14.84, 16.54, 
25.50. 27,34. 29.19. 31.06, 32,93. 
42.41. 44.32, 46 .23# 48.15. 50.05. 
59.73, 61.64, 63,60. 65,54. 67.49. 
77.23, 79.18, 81.14, 83.10. 85.06/ 


da t a ( hk ( n « 1 7 ) , 

1 4.37. 6.93, 

2 28.22, 31.25, 

3 57.43, 60.84. 

4 88.57, 92.08. 

5 120.32, 123.58, 


daU ( Hk ( n , 2 0 > . 

1 5.14. 7.40. 

2 31.51, 35,02. 

3 65.19. 69,14. 

4 Mil . 33. 105.40, 

5 138.24, 142.37. 


1 2.71. 

2 12.99. 

3 24.74, 

4 36.99. 

5 49.44. 


3.71. 4.74. 5.85 
14.25. 15.53. 16.82 
26.09, 27.44, 24.79 
34.34, 39.74, 41.12 


3-80. 4.44. 

14.34. 15.44, 


14.34, 15.71, 17.01, 14.46 
<*4 84. 24.26. 29.7 0, 31.13 
39 90. 4 1-36. 4 2.87. 44 . 37 
53.52, 55.07. 56.63, 58.20 


{Hklr>,25 > ,«M . 51 )/-2 • 394, -1 .60S. - 0 . 805. 0 .016 









DATA (HKCN.33). 

.MM , 51 )/-? . 376 . -1 . 582 

,-0.786 

. 0.015 

, ft . 8 4 2 

» 1^715 

939 

,-2.100 

.-1101 

, 0 . 023 

. 1.394 

, 3.015, 


I 2,64, 

3 .6? 

4.65, 5.73. 6.85, 

8 02. 

9.22, 

10.4*. 

11.73, 

55, 

14.08, 

U . 72 . 

19.47, 

22.31, 

25.23, 


2 1503, 

14.57 

15.75. 1 ' . 1 2 . 13.54, 

14.97, 

21.45, 

22.91. 

?4 .40. 

58 , 

40.80, 

44.06, 

47.36, 

50.69, 

54.05, 


3 25.91, 

2 7.45 

28.97, 30,51. 32.07, 

53.6 3. 

35.21 . 

34 .78. 

38.37. 

70 . 

7116. 

74.62, 

78.10, 

8 1.58, 

85.07 , 


4 39,96. 

41.55 

43.14, 4. 73. 46 .33, 

4 7.93, 

49. 52. 

51.11, 

52.7 1. 

1 1 , 

1 02 .43. 

104.16. 

109.69. 

113.23. 

116.77 , 


5 54,30. 

55.89 

57.47, 59.05, 60.63. 

62. 21 , 

63.78. 

65.35, 

66 .92/ 

01 . 

134 . 59 , 

138.17. 

141 .76, 

145.36, 

14 8 . 96/ 

C 

♦ • ♦ 













c 

c 

»** KFI02 














DATA ( HK t II . 34 ) , 

MM. 51 )/- 2. 772.-1 .9 28 

. -1 . 016 

, 0.021 

> 1.197 

. 2.496, 

842 

.-2.033 

.-1.104 

, 0.023 

. 1 . *30 

■ 3.138, 


1 3 - 90 . 

5.38. 

6.92. 8.51. 10.15, 

11.82. 

13.52, 

15.25, 

17.01, 

55. 

15-39. 

18.39. 

21.51. 

24.75, 

2 3.08, 


2 18 . 79 , 

20.58. 

22.40, 2'*.2J, 26,07 . 

27.93. 

24-89. 

31.68. 

3 3.57, 

24 , 

4 5.9 7 . 

49.73, 

53.54. 

57.39 , 

6 1.27 , 


3 35.47. 

37,37. 

59.28, 4!.|9, 43.11, 

45.03, 

46.96. 

48.69. 

5 0.82. 

09 , 

81-10. 

85.12. 

8 9. 16, 

93.20, 

9 7.26. 


4 52.7 5, 

54.69, 

54.62, 58.56. 60.51. 

62.45, 

6 4.39 . 

46.34. 

68.29. 

. 57 , 

117.64. 

121.77. 

1 25.87. 

129.99, 

134.11, 


5 70.24, 

72.19, 

74.14, U, . 1 0 . 78.05, 

60.01, 

81.97. 

83.94, 

8 5.90/ 

, 66 . 

154.82, 

1 58.99 , 

163.16, 

1 6 7.34. 

171. 52/ 

c 

*** 













c 

*** HNQ 3 














DATA (HKMI.351. 

71 - l . 51 IS-?. 83ft, -2.041 

,-1.139 

. 0 . 024 

, 1.442 

. 3.071, 








1 4,87. 

6.80, 

8.83, 10.94, 13. 1ft, 

15.31. 

17.56. 

19.83. 

22.13. 

. 393 

,, -2.487 

,-1.403 

. 0.030 

. 1 . A 98 

, 4.200, 


2 24.45, 

26.80. 

29.16, 31 .53, 33.92, 

36 .35. 

38.74. 

4 1.16. 

4 1 . 58 . 

.84, 

20.61 . 

24.58, 

28.69, 

32.93, 

37.30, 


3 46.02. 

48.45, 

50.90, 55.34, 55.79, 

58 . 23. 

60.68. 

63.13. 

65.51. 

.71. 

40.50. 

65-35. 

70.26. 

75.21 . 

80.19, 


4 68.04. 

70.49, 

72.94. 75.54. 77.65, 

60.30. 

82.76, 

85.21 . 

8 7 6 7 . 

.45. 

1 115. 54. 

110.67, 

115.81 . 

120.96,126.11, 


5 90.13. 

92.59, 

95.05, 97.52, 99.98. 

1 02.45 , 

1 04.92 . 

107 . 39. 

109.86/ 

.sz’. 

152.01 . 

157.21, 

162.41. 

167.62, 

1 7 2.8 4 , 

c 








.7 7 . 

199.03. 

204.29, 

209 . 56. 

214.84, 

220.13/ 

c 

» H202 














DATA ( HK ( 

II. 36 ), 

HM . 51 >/ - ? . 667 , « 1 . 845 

, -0 . 952 

, 0.019 

, 1 . 08 0 

, 2.254, 








1 3.47, 

4.79, 

6.16, 7.53, 9.05. 

10.56. 

12 .11. 

13.69. 

15.31, 

. 387 

,-l . 593 

,-0 . 796 

, 0.015 

. 0.864 

, 1 . 760 , 


2 16.95, 

18.62, 

20.32. 22.04. 23.77, 

25.53. 

27.31, 

29.10. 

30.90. 

.85. 

6.98, 

8.14. 

9.32, 

10.53, 

11.75, 


5 32.72, 

34.55, 

36.40, 38.25, 40.11, 

41.98. 

43.56, 

45,75. 

47.64, 

.82, 

13.11. 

19.42, 

20.74, 

22.07 , 

23.40, 


4 49.54. 

M . 4 5 . 

53.36, 55.28, 37.21, 

59.14. 

61 . 07 . 

63.01, 

64.96, 

. 79, 

30.15. 

31.51. 

32.88. 

34.25, 

35.62, 


5 66.91, 

68 .86, 

70.63. 72.79, 74.76, 

76.74. 

78.72. 

60.70, 

82 69/ 

.12. 

42.50. 

6 3.89. 

45.27 . 

46.66, 

48.05. 

c 

• » + 







.64, 

55.04. 

56.44, 

57 .85. 

59.26, 

6 0 .67/ 

c 

*♦» iica 














DATA < MKC IT. 37 ) . 

II* 1 . 51 )/~2. 44 3,-1 .691 , 

, -0.890 

. 0.016 

, 1.024, 

, 2.111, 








1 3.26, 

4.47, 

5.73. 7.02, 8.34, 

9.69, 

11.05. 

12.42. 

13.81, 

. 392 

:.-i .59* 

, -0.800 

, 0.015 

, 0.S78 

, 1.798, 


2 15.21. 

16.63, 

18.05, 19.48. 20.92. 

22.36. 

23.81 . 

25.26 . 

24.72, 

. 94 . 

7.09, 

8.24, 

9.42. 

10.62. 

11.85, 


3 28.18. 

29.64, 

31.11. 32.57. 34.04. 

35.51. 

36 . 98 . 

38 .45. 

39 .92, 

.95, 

18.27. 

19.60, 

20.95, 

22.31. 

23 .69. 


4 41.40, 

4 2.8 7, 

44.34, 45.81, 47.29, 

48 . 74 . 

50.24, 

51.71. 

53.19, 

. SI. 

30.74. 

32.13. 

33.6 3, 

35.08. 

36 . 54, 


5 54.66, 

56.14, 

57.62, 59.10, 60.57, 

62.06, 

*3.54, 

65 . 02. 

66.50/ 

. 4 J. 

43.92, 

45.40. 

46 .89 . 

48 .39, 

49.88. 

c 

♦ * » 







.87, 

57.37, 

58.86, 

60.56, 

61.86. 

63 . 56/ 

c 

♦ ** IIH 













c 

DATA (HKCN.38). 

11= 1 . 51 )/- 2 . 0 56 . -1 . 334 „ 

, -0 . 685, 

, 0.013, 

, 0.710, 

, 1 .409 . 








1 2.11, 

2.82, 

3.54, 4.26. 5,01. 

5.76. 

6.53. 

7.31, 

6 .10. 

46 0 

.-1.645 

.-0.839 

. 0.016 

, 0.944 

, 1.949, 


2 8.90. 

9.71, 

10.53. 11.36. 12.20, 

1 3 05. 

1 3 90. 

14.76. 

15.63, 

57 . 

7 .82, 

9.10, 

10.34, 

11.70, 

13.02. 


3 16.51, 

17.59, 

18.28. 19,1.8. 20.08, 

20.99, 

21.91. 

22.83, 

23 75. 

46 . 

19 . A2. 

21 .21 . 

22.60. 

24,01. 

25.42. 


4 24.69. 

25.63. 

26.57 . 27.52, 23.4.8. 

29 . 44. 

30.41 , 

31.38, 

32.36. 

13. 

32.53. 

34.01, 

35.49. 

34 . 75. 

38 . 42 , 


5 55.34. 

34 . 53, 

35.33, 56.33. 37.33, 

55.34, 

39.36, 

40.38 . 

4 l . 40/ 

37 . 

4 5 . * 7 . 

47.39. 

48.91. 

50.43, 

5117. 

c 








20 . 

59.79, 

61.58. 

62. 9 7 . 

64.58. 

66 .20/ 

c 

IIH2 













: c 

DATA ( HK ( II , 39 ) • 1 

HM . 51 >/- 2 . 37 5, - l . 387 . 

-0 .788. 

, 0.015, 

0.834, 

, 1.679. 


4 '.s! 5:95. 7!53. 9.24, 11.04, 12.94. 14.91, 16.96, 

2i:24 , 2 3.50, 25.80, 28.14, 30.52. 32.94. 35.40, 37 .90, 


40 42. 42.98. 45.56. 
4419. 66.9 2, 49.67, 
89.30, 92.16, 15.04. 


48 14, 30.74, 53.43. 56.09. 58.77. 61.47, 
72.43. 75.20. 77.99. 80.80. *3.62. 86.45. 
97. 9 <•. 100. 84. 103. 77. 106. 71. 109. 66, 112. 6 3/ 


.26). a*1. 51 )/-2. 391,-1 696.-0,138. 0.019, 1.151 
5.I8. 4.66. 8.20, 9.79. 1 1.44. 13. !<.. 14,85. 


4ft. 2ft. 23. 22 08, 23.94 
42. 37 .61 . 39.62. 41 .64 

92. 5 5.99. 58.07, 60. 16 


25,85. 27.77. 29 


51 .66 . 33 . 6 3. 


<• 5.70. 47.74, 49 79. 51 .85. 


7495 ! 77.09, 79.23, 81.38. 83.54. 85.70, 87.87, 90.05/ 


DATA I HA< M. 40 ) . M* 1 . 5 1J/-2 . 790 


, - 1 . 594 , - 0 . 
818 . 9 . 

20 . 74 . 22 . 
34 12 . 35 . 
47 . 78 , 49 
61 . 65 . 63 . 
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APPENDIX C 

INPUT FILE FOR SIMPLIFIED MECHANISM WITH FIVE CHEMICAL REACTIONS 


T3AAA KOI I 887 13X1X20 


irest 

0 

i p os t 

1 

n x 

13 

n y 

1 

n 2 

20 

1 u/a 1 1 

1 

nc hop 

0 

1 pr 

0 

j sec tr 

1 

irez 

0 

nc f i lm 

200 

nc tap8 

9999 

nc 1 os t 

4000 

c a f i 1 m 

5. 0 

c a f i n 

360. 0 

cadump 

-90. 0 

dcadmp 

10. 0 

a ng mom 

1. 0 

cy 1 

1. 0 

dy 

0. 0 

p g s s ui 

1. 0 

samp 1 

0. 0 

dti 1 

. 04167«-4 

dtmxca 

9. 97e+9 

d tma x 

9. 99e+9 

1 1 i md 

0. 0 

twf i 1 m 

9. 99e +9 

tw f i n 

9. 99e+9 

f c h sp 

0. 25 

stroke 

0. 0 

squi s h 

100. 0 

rpm 

0. 0 

a t d c 

0. 0 

c onr od 

16. 269 

offset 

0. 0 

su/i r 1 

0. 0 

sun pro 

1. 01e-10 

t h s ec t 

0. 5 

ep sy 

1. 0e-3 

ep s v 

1. 0e-3 

epsp 

1. 0e-4 

e p s t 

1. 0e-3 

ep s k 

1. Oe— 3 

e p s e 

1. Oe— 3 

9 ^ 

0. 0 

gy 

0. 0 

9* 

0. 0 

t c y 1 uil 

400. 0 

thead 

400. 0 

t p i s tn 

400. 0 

tva 1 ve 

400. 0 

TEMPI 

400. 0 

pardon 

0. 0 

aO 

0. 1 

bO 

1 . 0 

an c 4 

0. 05 

a d i a 

0. 0 

anuO 

0. 0 


visrat-. 6666*667 
tcut 800. 0 
tcu te 1200. O 


• p ichni 


omg c hm 
TKEI 


0. :*o 

1 . o 

2. 2SOE+4 


2-D 


COMBUSTOR a 


TKESW 

1. 0 

s g s 1 

0. 0 

un i sea 1 

0. 0 

a irmu 1 

i. 457e-5 

a i rmu2 

i 10. 0 

a i r 1 a 1 

252. 0 

a i r 1 a2 

200. 0 

e x p d i f 

0. 6 

pr 1 

0. 74 

r pr 

1.11 

rprq 

1. 0 

r p r e 

0. 769231 

r sc 

1.11 

x i g n i t 

2. 0e+3 

T1IGN 

0. 0 

TDIGN 

9. 99E+9 

c a 1 i gn 

-27. 0 

cad i gn 

9. 6 

i i gnl 1 

1 

i i g nr 1 

1 

j i gn f 1 

1 

j i gnd 1 

1 

k i g nb 1 

2 

k i gn 1 1 

n 

i i g n 1 2 

1 

i i gnr2 

1 

jig n f 2 

1 

j i g n d 2 

1 

k i gnb2 

"3 

k i g n t2 

n 

k wi k e q 

0 

numn o z 

1 

numve 1 

1 

i n j d i s t 

1 

k o 1 i d e 

1 

1 1 i n j 

0. 0 

t d i n j 

9. 99e + 9 

ca 1 i n j 

-52. 0 

cad i n j 

12. 672 

t spinas 

9. 208 

pulse 

0. 0 

tnparc 

4000. 0 

r h op 

0. 7436 

t p i 

350. 0 

t ur b 

1 . 0 

breakup 

1. 0 

e vap p 

1. 0 

d r no i 

0. 0 

d z no z 

0. 1 

d t lino z 

0. 0 

t i 1 1 x y 

0. 0 

t i 1 1 x z 

0. 0 

cone 

62. 5 

d c one 

12. 5 

ano z 

1. 0 

smr 

5. 00e-4 

amp 0 

0. 0 

4000. 0 

np o 

14 

nun i f 

0 

1 1 

0. 0 

2 1 

. 195 

3 1 

. 391 

4 X 

. 5P6 

5 1 

. 702 


0 . 0 
0. 0 
0. 0 
o. o 

o. o 
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APPENDIX C - Continued. 


b 1 . 977 0. 0 

7 1 1. 172 0. 0 

Q 11. 360 0. 0 

9 11. 563 0. 0 

10 1 I- 75B 0. 0 

11 11. 954 0. 0 

12 1 2- 147 0. 0 

13 1 2 . 345 0 . 0 

14 1 2. 540 0. 0 


nho 0 

s qua re 0.0 

r c or nr 0. 0 

ns t r t 0 


icont 1 1101 1 1 101 1 1 1000000001 1000 


mirror 


n v zone 

n 

c_ 







0. 0 

-5. 0 

1. OelO 0 

0 0 

0 

0 1 0 



-5. 0 

-45. 0 

36. 0 0 

1 1 

1 

1 1 0 



n vvve c 

4 







n vp ve c 

0 







0. 0 

-45. 0 

15. 0 

0 

1 

0 



3. 0 

-45. 0 

15. 0 

0 

3 

0 



0. 0 

-5. 0 

1 . OelO 

0 

0 

6 



0. 0 

-5. 0 

1. OelO 

0 

0 

13 



n vc o n t 

n 







0. 0 

-45. 0 

1 5. 0 

0 

3 

0 



0. 0 

-5. 0 

1. OelO 

0 

0 

6 



n s p 

12 







c3h0 

r h o 1 

0. 0 






o 2 

rho 2 

3. 1096e-4 

mw 2 


32. 000 

h t f 2 

0. 0 

n 2 

rho3 

1. O0O9e-3 

mw3 


28. 016 

ht f3 

0. 0 

c o 2 

r h o 4 

1 . 4 007e -5 

mw4 


44. 01 1 

h t f4 

-93. 965 

h 2 o 

r h o 5 

7. 0030e -6 

mui5 


18. 016 

ht f5 

-57. 103 

h 

r h 06 

0. 0 

rmjj 6 


1. 000 

h t f 6 

51. 631 

h 2 

r h o 7 

0. 0 

mu7 


2 . 016 

htf7 

0. 0 

o 

rh oQ 

0. 0 

niu 0 


16. 000 

h t f 8 

50. 909 

p 

r h o 9 

0. 0 

mw9 


14. 008 

h t f9 

1 12. 520 

o h 

r h o 1 0 

0. 0 

mui 1 0 

17. 008 

h tf 10 

9. 239 

c o 

r h o 1 1 

0. 0 

mw 1 1 

20. 01 1 

h t f 1 1 

-27. 200 

n o 

r h o 1 2 

0. 0 

muj 1 i 

2 

30. 008 

ht f 12 

21. 456 


r t ou t 0.0 

topout 1.0 

b o t i n 1.0 

d i s tamb O. O 
pamb 1.6229e+6 
t k eamb 0.10 

sc lamb 0. 273450 
spdambl 0.0 

s p d amb 2 3. 1096e-4 
spdamb3 1. 0809e— 3 
s p d amb 4 1. 4007e— 5 

s p damb 5 7. 0030e-6 
spdamb6 0. 0 
spdamb7 0. 0 
s p d amb 8 0. 0 

s p d amb 9 0.0 

s p d a m 1 0 0. O 
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TABLE I. - ENGINE PARAMETERS AND TEST CONDITIONS 


Engine stroke, cm 9.55 

Compression ratio . . . . 6.54 

Volume of piston bowl, cm^ 58.8 

Squish clearance, cm 0.1819 

Connecting rod length, cm 16.269 

rpm 1600 

Initial engine air temperature, K 400 

Initial engine pressure, atm 1 

Fuel ..... isoctane 

Starting crankangle of injection, deg 52 BTDC 

Injection duration, deg 12.672 

Injection angle, deg 60 

Starting crankangle of ignition, deg 27 BTDC 

Ignition duration, deg 9.6 
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Mr...... ... V X \ \ \ \ 1 i 



U UAV - 7.89828E+02, V MAV = 2.94169E+03, 309 PARTICLES IN THE SYSTEM 

MAX MAX 

W MAV - 9.76695E+02 

MAX 

(a) GAS VELOCITY AT 29.05° BTDC. (b> FUEL SPRAY AT 29.95° BT1X . 


M = 8.84606E+02, V M ,„ = 3.38367E+03. 122 PARTICLES IN THE SYSTEM 

MAX 

W MAV = 1.245S5E+03 

MAX 

(C) GAS VELOCITY AT 11.37° BTDC. (d) FUEL SPRAY AT 11.37° BTDC. 



L = 6.57971E+02, H = 2.9607^E+03, c 8 h 18 ACR0SS J = 1 PLANE ' L ~ 1.17529E-02, 

DQ - 2.30277E+02 H = 1.05776E-01, 11.37°, DQ = 1.17529E-02 

<e> TEMPERATURE AT 11.37° BTDC. (f) FUEL VAPOR AT 11.37° BTDC. 

FIGURE 1. - RESULTS OF DISC GASOLINE-FUELED PISTON ENGINE FOR CASE 1. 
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VELOCITY ACROSS J = 1 PLANE, -29.17°, 472 PARTICLES IN THE SYSTEM 

U max = 7.86088E+02, V^ = 2.94103E+03, 

W max = 9-79231E+02 

(a) GAS VELOCITY AT 29.17° BTDC. (b) FUEL SPRAY AT 29.17° BTDC. 




L = 4.83900E+02, H * 6.50552E+02, 
DO = 1.66652E+01 

(c) TEMPERATURE AT 29.17° BTDC. 


U hax = 1-02988E+O3, = 3 . 34960E+03 , 

W HAX = 1 ■ 08345E+03 

(d) GAS VELOCITY AT 13° BTDC. 





432 PARTICLES IN THE SYSTEM 

(e) FUEL SPRAY AT 13.89° BTDC. 

FIGURE 2. - RESULTS OF DISC GASOLINE 


L = 6.82538E+02, H = 2.92452E+03, 
DQ = 2.24198E+02 

(f> TEMPERATURE AT 13.89° BTDC. 
-FUELED PISTON ENGINE FOR CASE 2. 


L: L V * -uji i i* • i m I ‘ /cCi 51 ft 5 

OF POOR QUALITY 
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L = 4.83019E-07, H = 1.91309E-02, 
DQ = 1.91313E-03 
(g) FUEL VAPOR AT 13.89° BTDC. 


L = 0 . OOOOOE+OO , H = 2 . 1 3202E-K)0 , 
DQ = 2.13202E-01 


(h) C0 2 AT 13.89° BTDC. 


FIGURE 2. - CONCLUDED. 
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II . / 

VELOCITY ACROSS J = 1 PLANE, -1.82° 
U MAX = 7.02298E+02, V MAX = 1.69386E+03, 

W MAV = 1 .2761 1E+03 

MAX 

(a) GAS VELOCITY AT 1.82° BTDC. 



123 PARTICLES IN THE SYSTEM, -1.82° 
(b) FUEL SPRAY AT 1.82° BTDC. 



.TEMP ACROSS J = 1 PLANE, L = 9 . 19304E+02, 
H = 2.98323E+03, -1.82°, DQ = 2.29325E+02 

(C) TEMPERATURE AT 1.82° BTDC. 



c ? h 8 ACROSS J = 1 PLANE, L = 7.23420E-03, 
H = 6.51078E-02, -1.82°, DQ = 7.23420E-03 

(d) FUEL VAPOR AT 1.82° BTDC. 


FIGURE 3. - RESULTS Of DISC GASOLINE-FUELED PISTON ENGINE FOR CASE 3. 
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VELOCITY ACROSS J = 1 PLANE. -0.16°, 28 PARTICLES IN THE SYSTEM 

U HAX = 1 • 20566E+03 , V^ = 3.59255E+03. 

H «ax = 1 • 16769E+03 


(a) GAS VELOCITY AT 0.16° BTDC. (b) FUEL SPRAY AT 0.16° BTDC. 



■ TEMP ACROSS J - 1 PLANE, L - 9.71871E+02. Cj g ACROSS J = 1 PLANE. L = 9.09029E-03, 

H = 2.62905E+03, -0.16°. DO - 1.65718E+02 H = 8.18126E-02, -0.16°. DO = 9.09029E-03 

(C) TEMPERATURE AT 0.16° BTDC. (d) FUEL VAPOR AT 0.16° BTDC. 

FIGURE 9. - RESULTS OF DISC GASOLINE-FUELED PISTON ENGINE FOR CASE 9. 
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VELOCITY ACROSS J = 1 PLANE, VELOCITY ACROSS J = 1 PLANE, T = 4.15503E-03, 

T = 8.26447E-03, CYCLE 1200, CYCLE 1200, U MAX = 5.59223E+04, V^ = 

U MAX - 1 . 20066E+04 , V MAX - 1.21813E+04, 2 .S3656E+04, W MAX - 5.33020E+04 


Vx = 1 .73781E+04 



TEMP ACROSS J = 1 PLANE, L = 1 .07907E+03, TEMP ACROSS J 1 PLANE, L 5.29182E+02, 

H = 1 58762E+03 , T = 8.26447E-03, CYCLE H = 1.48754E+03, T = 4.15503-03, CYCLE 
1200, MIN = 1 .01551E+03, MAX = 1.65119E+03, 1200, MIN = 4.09387E+02, MAX - 1.60733E+03, 

DQ = 6 . 35683E+01 00 = 1 - 19794E+02 

(a) OCTANE COMBUSTION. <b) PROPANE COMBUSTION. 

FIGURE 5. - SIMULATION OF COMBUSTION IN TWO-DIMENSIONAL COMBUSTOR. 
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